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1.0 SUMMARY 


Program L221 (TEV156) is structured as four overlays, one main and three primary. 
Input into the program is made via cards and magnetic files (tapes or disks). Output 
from the program consists of printed results and magnetic files containing vectors 
suitable for plotting. 

Although L221 serves as a module of the DYLOFLEX system, it can be operated as a 
standalone program. Subroutines used by L221 include routines embedded in the 
program code, routines obtained from the standard FORTRAN library, and routines 
obtained from the DYLOFLEX alternate library. 


2.0 INTRODUCTION 


Program L221 (TEV156) was developed for use as either a standalone program or as a 
module of a program system called DYLOFLEX (see fig. 1), developed for NASA under 
contract NAS1-13918 (ref. 1). Because of the DYLOFLEX contract requirements 
developed in reference 2, a program was needed to calculate power spectral 
density (PSD) gust load parameters for equations of motion and load equations that 
have frequency -dependent coefficients. An existing program 1 that calculated PSD gust 
load parameters for equations with constant coefficients was modified according to the 
DYLOFLEX specifications 2 to also use nonconstant coefficients matrices for the 
equations of motion and load equations. 

The objective of this volume is to aid those persons who will maintain and/or modify the 
program in the future. To meet this objective, the following items are defined in some 
detail: 

• Program design and structure 

• Overlay purpose and description 

• Input, output, and internal data base descriptions 

• Test cases and procedures 

1 Clemmons, R. E.: A Power Spectral Digital Computer Program to Determine Dynamic Loads Due 

to Random Gusts-PSDSYS (TEV156) - Users Guide. BCS-G025, June 1973. 

2 

Clemmons, R. E.: Programming Specifications for Modules of the Dynamic Loads System to 
Interface with FLEXSTAB . NASA contract NAS1-13918. BCS-G0701, September 1975. 

(Internal Documents.) 




Figure 1.-DYL0FLEX Flow Chart 










3.0 PROGRAM DESIGN AND STRUCTURE 


The program is structured as a system of four overlays (fig. 2). 

Main overlay (L221,0,0) L221vc 
Primary overlay (L22 1,1,0) FINDRMS 
Primary overlay (L22 1,2,0) SORTQLS 
Primary overlay (L22 1,3,0) PLOTQLS 

The main overlay L221 reads cards that direct the execution of the primary overlays. It 
also aids communication between the primary overlays via labelled common blocks. 

The 1,0 primary overlay FINDRMS performs all of the analytical calculations of the 
program. Optionally, its results are written onto magnetic files (tape or disk) to be 
processed by the two other primary overlays. 

The 2,0 primary overlay SORTQLS sorts the FINDRMS results (responses - Q’s, loads, 
output spectrum, and A values) and prepares a magnetic file (tape or disk) of data 
vectors suitable for plotting. 

The 3,0 primary overlay PLOTQLS sorts the FINDRMS results exactly the same as 
SORTQLS, but writes a magnetic file of plot data and plotting instructions geared 
specifically to the COMp80 plotter. 

The input and output data for L221 (TEV156) is displayed in figure 3. Each overlay is 
discussed in detail in succeeding sections. Included for each overlay are: 

1. The overlay’s purpose 

2. The overlay’s analytical steps 

3. The input/output devices used 

4. A macro flow chart 

5. Table of subroutines called 

(Note: All subroutines have only one entry point) 

Special symbols are used on the tables of subroutines called to indicate routines that 
are loaded from the operating system library and from the DYLOFLEX alternate 
subroutine library (DYLIB). All other subroutines are local to L221 (TEV156). 

This document does not contain a description of each local subroutine. Please see the 
program listing, where each routine contains "internal documentation,” comments 
describing its purpose, operation, input, output, modification history, etc. 

The storage and handling of data is discussed in section 3.5, Data Bases. All magnetic 
file (tape or disk) communication outside the program is done via READTP/ WRTETP 1 
To achieve the most efficient use of core, a scheme to allow variably dimensioned arrays 
(dynamic storage allocation) was used. Section 3.5 also displays the variably 
dimensioned arrays (matrices) and describes the subroutines that keep track of them. 

1 Clemmons, R. E.: Programming Specifications for Modules of the Dynamic Loads System to 
Interface with FLEXSTAB. NASA contract NASl-13918, BCS-G0701, September 1975. 

(Internal Document.) 
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Control the execution 
of L221 's primary 
overlays 



Calculate the 
frequency responses, 
RMS loads, and 
number of zero 
crossings 


Sort frequency 
responses and 
prepare a magnetic 
file of data vectors 
suitable for 
plotting 


Sort frequency 
responses and 
prepare a magnetic 
file of plot instructions 
for the COMp80 plotter 


Figure 2.— Overlay Structure of L221 (TEV156) 



IPLTPE IRTAPE NEWTPE TAPE99 


Figure 3.— External Input/Output of L221 ( TEV156 ) 
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3.1 OVERLAY (L221,0,0) - L221vc 


The main overlay of L221 (TEV156) is L221vc, where v is a letter indicating the 
program version, and c is an integer number indicating the correction that applies to 
the v version. 

Purpose of L221vc 

L221vc performs certain bookkeeping tasks, directs the execution of the primary 
overlays, and aids communication between primary overlays via labeled common blocks. 

Analytical Steps of L221vc 
L221vc performs its task in the following Bteps: 

1. The subroutine PRGBEG is called to place the program header on the printed 
output. 

2. A data card is read. It must begin with $FREQ to assure that the card input file is 
correctly positioned. If it does not contain $FREQ, execution is terminated. 

3. A program directive card is read, printed, interpreted, and acted upon according to 
the following logic: 

a. If the keyword is $TITLE jump to step 3 again. 

b. If the keyword is $CHECKOUT, set the checkout switch and jump to step 3 
again 

c. If the keyword is $FIND, jump to step 4. 

d. If the keyword is $SORT, jump to step 5. 

e. If the keyword is $PLOT, jump to step 6. 

f. If the keyword is $QUIT jump to step 7. 

4. Overlay (L221,l,0) is called. When it is finished program control returns to step 3. 

5. Overlay (L221,2,0) is called. Wl>en it is finished program control returns to step 3. 

6. Overlay (L221,3,0) is called. When it is finished program control returns to step 3. 

7. Subroutine PRGEND is called. This subroutine places the program termination 

message on the printed output. 

If fatal errors are discovered, the program prints a diagnostic and jumps to step 7. The 
macro flow chart of this overlay is shown in figure 4. The subroutines called are 
displayed in table 1. 


I/O Devices of L221vc 

L221vc reads program directive cards and writes them on the printed output file. All 
other I/O accomplished by L221 (TEV156) is done in the primary overlays. 


5 




Program 
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► print an 

input card image 
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Figure 4.— Macro Flow Chart of Overlay (L221,0,0)—L221vc 
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Read and 
print 

input card 
image 


Obtains 
the control 
card keyword 



Figure 4. —(Continued) 
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Figure 4.— (Concluded) 







Table 1. —Routines Called by L221vc 

OVERLAY (L 221, 0,0) 

PROGRAM L221VC 


CHKWRD 
FINDRMS 
IRDCRD* 
PLOTQLS 
PRGBEG+ 
PRGEND + 
REQ FL + 
SORTQLS 


OVERLAY (L2 21 ,1,0) 
OVERLAY (L221 ,3,0) 


OVERLAY (L221 ,2,0) 


+ indicates a routine on the DYLOFLEX alternate 
subroutine library DYLIB. 

All others are local to L221 (TEV1 56) . 



3.2 OVERLAY (L221,l,0) - FINDRMS 


Purpose of FINDRMS 

The first primary overlay of L221 (TEV156) is FINDRMS. FINDRMS performs all of the 
analytical tasks of L221. Given an input power spectrum </>(Q) over a range of 
frequencies (Dj, £^2, ...» flNFREQ)> matrices of equations of motion and load equations, 
and excitation forces, FINDRMS will find: 

{Q} The generalized coordinates at each frequency 

{SUMC} The load transfer function at each frequency 

and integrate over the range of frequencies to find: 

{A} The RMS load values 

{N 0 } The number of zero crossings per unit length 

Analytical Steps of FINDRMS 

The FINDRMS operations are divided into the following steps: 

1. Read card input of options and constants. 

2. Read card or tape input of equations of motion ([Mj], [M2 ], [Mg] and [S]], 

S2 ], [S 6 ] for feedback) plus the excitation forces ({ C2 } and {C3} or {f/} and 

</>]). This includes {FREQM} frequencies if NKVAL > 0. 

3. For each of the NFREQ frequencies, solve for { Q} , the generalized coordinates, 
save {Q} on the scratch file, and optionally print { Q} . 

4. (Optional) - Modify the equations of motion for a static elastic solution and repeat 
step 3. Chosen columns of the matrices [M2], [M3], [M5] and [Mg] are set equal 
to zero for the static elastic solution. 

5. Read card_or tape input of load equations ([Mj], [M2],_..., [Mg] and [Si], [S 2 ], 
[S 3 ], ..., [Sg] for feedback) plus the matrices {C2} and {C3} or [</>]. 

6. For each frequency read {Q} (from step 3), calculate {SUMC}, the input spectrum 
</>(£!), and the output spe ctrum SPEC, optionally write data on files IPLTPE and 
IRTAPE, optionally print {SUMC}, and keep a running integration over the frequen- 
cies of {A} and {N 0 }. 

7. Print {A] and {N 0 } and optionally punch them on cards. 

8. (Optional) - Check the correlation between the different load responses requested 
by card input data. 

9. (Optional) - Modify the load equations for a static elastic solution (set chosen 
columns of [M2], [M3], [M5] and [Mg] to zero) and repeat steps 6, 7, and 8. 

10. (Optional) - For additional sets of load equations repeat steps 5 through 9. 

Figure 5 contains a macro flow chart of FINDRMS. The subroutines called by FINDRMS 
are displayed in table 2. 

I/O Devices of FINDRMS 

The possible I/O devices for FINDRMS are shown in figure 6. For a complete description 
of the input data cards and magnetic files used, see sections 6.3 and 6.4 in volume I of 
this document. 

The two scratch files, SCRATCH and SCRAT2, are described in section 3.5.3. 
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of F1NDRMS 


Initialize 
- theVARDIM 
array storage 

* Read Block I 
card input 
data, check for 
« fatal errors, 
and interpolate 
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Initialize the 
output and 
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Figure 5. —Macro Flow Chart of Overlay IL221 , 1,0)—FINDRMS 
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{ Read and interpolate for 
frequency-dependent 
matrices and check 
for any fatal errors 


Calculate {C 3 | 
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degrees of freedom 
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Check for end-of- 
frequency loop 


Figure 5.— (Continued) 
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Is the static- 
elastic solution 
required? 
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been done? 


r Modify the 
equations of 
- motions for a 
static-elastic 
k solution 


Read Block 1 1 1 
card input 
data and 
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end of input data 


Figure 5.— (Continued) 
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Figure 5 . -(Continued) 


15 









Figure 5.— (Continued) 
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Figure 5 .— (Continued) 
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Release the 
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f 
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termination 
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Figure 5.— (Concluded) 
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Table 2.— Routines Called by FI NORMS 

OVERLAY (L221, 1 , 0) 

PROGRAM FINDRMS 


DELETR+ 

FETADD+ 

FETDEL + 

FRCALC 

CONDAT INITIR+ 

LOCF* 

NAMFIL+ 

READTP+ 

REQFL+ 

* 

f ENDFIL* 

SINTRP J FETADD+ 

FETDEL+ 

* indicates a routine in the FORTRAN subroutine library. 

* indicates a routine in the DYLOFLEX alternate subroutine library DYLI B . 

All others are local to L221 (TEV156). 
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Table 2. -(Continued) 


OVERLAY (L2 2 1 , 1,0) 
PROGRAM FINDRMS 


CORREL 

FETADD+ 

FETDEL+ 


GCMAT 


DELETR+ 
FETADD+ 
INITIR+ 
< TRDCRD + 
LOCF* 

MATRED 


REQFL+ 
READTP ♦ 


EOF* 

INITIR+ 

< 

PRINTM+ 

PRNTCM+ 

READTP+ 


GRDPEN 



Table 2.—( Continued) 


OVERLAY (L22 1 # 1 , 0) 
PROGRAM FINDRMS 


LOMAT 



DELETR* 
EOF* 
FETADD * 
FETDEL+ 
FSF + 
IRDCRD+ 

INITIR+ 

LOCF* 

MATH ED 


NAMFIL+ 

REQFL+ 

R EADTP + 


EOF* 

INITIR+ 

PRINTM* 
PR I NTCM + 
REA DTP* 
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Table 2.— (Concluded) 


OVERLAY (L221 r 1 ,0) 
PROGRAM FINDRMS 


MATINT 

OUTAP 

RETURN F + 
SECOND* 
SAV EQ 

SOLVEL 

SOLVEQ 

SPCTRM 
STARTR+ 
STATE L 


EOF* 

INITIR+ 

< MATRED < PRINTM+ 

L PRNTCM+ 

READTP+ 

< ENDFIL* 

WRTETP+ 


2 ENDFIL* 
FSF + 

J FEDBAK 
SUMATS 
"CGLESM+ 
1 FEDBAK 
SUMATS 
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i 1 1 mm min 


3.3 OVERLAY (L221,2,0) - SORTQLS 
Purpose of SORTQLS 


SORTQLS is called to process data written on the magnetic file IPLTPE by FINDRMS, 
and write a magnetic file, NEWTPE, which contains pairs of matrices, independent and 
dependent variables, to be plotted by a subsequent program. 

IPLTPE contains one file of data for each FINDRMS solution. Each file begins with the 
frequency array {FREQ} (a 1 x NFREQ matrix), which is followed by NFREQ groups of 
three or four of the following arrays ({A} will be present only if IPLRMS = 1 on FINDRMS 


card set 4.3). 



o {QMAGr} 

1 x NDOF 

The magnitude of Q 

o {QMAG l } 

1 x NLD 

The magnitude of SUMC 

o {SPEC} 

1 x NLD 

The output spectrum 

o {A} 

1 x NLD 

The RMS loads, through 
frequency i 


If the static elastic solution was requested in FINDRMS, every second file on IPLTPE is 
from a static elastic solution. 


Analytical Steps of SORTQLS 

The SORTQLS operations are divided into the following logic: 

1. Read a set of data cards defining the SORTQLS options^the FINDRMS load set to 
be processed, and the item (QMAGr, QMAGl, SPEC, or A) to be sorted. 

2. Position IPLTPE at the beginning of the FINDRMS solution to be processed. 

3. Read and store the frequency array, {FREQ}. 

4. Read and store the NFREQ matrices of the desired item (QMAGr, QMAGl, SPEC, 
or A) as the rows of the NFREQ x NR matrix [SOLl] shown below. 

SOLI 


NFREQ < 


Row 1 
Row 2 


where NR = NDOF for item = QMAGr 
and NR=NLD for item QMAGl, 

SPEC, or A. 




L Row NFREQ J 

Vi y ^ 

NR 
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5. (Optional) - When pairing of standard and static elastic solutions is requested, step 
4 will be repeated for the same item of the next solution, which is the static elastic 
solution in the subsequent file. The new solution is stored in the array [SOL2], 
which is identical to [SOLl] in size. 

6. (Optional) - Scale the frequencies and/or [SOLl] (and [SOL2]) as requested. 

7. Write on NEWTPE the frequency array {FREQ} and the column of [SOLl] 
corresponding to the chosen degree of freedom (for QMAGr) or load (for QMAGl, 
or SPEC, or A). 

8. (Optional) - Repeat step 7 for the array [SOL2] if pairing was requested. 

9. Repeat steps 7 and 8 for all generalized coordinates, loads, spectra, or A values 
requested by the current data set. 

10. Repeat steps 1 through 9 for as many data sets as needed, but terminate SORTQLS 
execution after reading the $END card (see card set 61). 

Figure 7 contains a macro flow chart of SORTQLS. The subroutines called by SORTQLS 
are displayed in table 3. 
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Initialization 
of SORTQLS 


Read a set of 
card input 
data and 
check for fatal 
errors 


Figure 7. —Macro Flow Chart of Overlay (L221 ,2,0)— SORTQLS 
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to sort 


Figure 7.— (Continued) 
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Figure 7. —( Concluded ) 
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Table 3.-Routines Called by SORTOLS 

OVERLAY (L221 ,2 ,0) 

PROGRAM SORTQLS 


DELETR+ 

ENDFIL* 

FETADD+ 

FETDEL+ 

TNITIR+ 


KARDIN 


CHKWRD 
< EOF* 
IRDCRD+ 
NAMFIL+ 

V- 


LOCF* 

REDSOL 

REQFL+ 

SECOND* 

STARTR+ 

WRTETP+ 


<|rEADTP* 


* indicates a routine in the FORTRAN subroutine library. 

♦ indicates a routine in the DYLOFLEX alternate subroutine 
library. 

All others are local to L221 (TEV156). 
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I/O Devices of SORTQLS 


The possible I/O devices for SORTQLS are shown in figure 8. For a complete description 
of the input data cards and magnetic files, see sections 6.3 and 6.4 in volume I of this 
document. 



NEWTPE 


Figure 8.— Overlay (1-221,2,0) SORTQLS Input/Output 
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3.4 OVERLAY (L221,3,0) - PLOTQLS 


Purpose of PLOTQLS 

The Overlay (L221,3,0), PLOTQLS, is called to process the data written on file IPLTPE 
by FINDRMS and write an output file TAPE99, which contains the plotting instructions 
for the COMp80 plotter. 

Analytical Steps of PLOTQLS 

PLOTQLS sorts the independent and dependent variables in exactly the same manner 
as program SORTQLS (see section 3.3). Steps 7 and 8 are replaced by calls to NPS 
subroutines, 1 which generate the plotting instructions. 

PLOTQLS has additional data instructions not applicable to SORTQLS; the 
establishment of grid limits and plot labeling information are two examples. 

Figure 9 contains a macro flow chart of PLOTQLS. The subroutines called by PLOTQLS 
are displayed in table 4. 


Note: The automatic plotting overlay PLOTQLS in this program requires subroutines 
that are proprietary to The Boeing Company. 


1 Numerical Plotting System - Users Manual . BCS-G0509, March 1976. 


(Internal Document.) 
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Initialization 
of PLOTQLS 


Check for 
previous calls 
to PLOTQLS 


Unblock 

TAPE99 


Initialize 
NPS package 


Figure 9.— Macro Flow Chart of Overlay (L221 ,3,0)— PLOTQLS 
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Figure 9. —( Continued ) 
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Figure 9. —( Continued ) 
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Figure 9 .— (Continued) 
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Figure 9. —( Concluded ) 
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Table 4,— Routines Called by PLOTQLS 


OVERLAY (L22 1,3,0) 
PROGRAM PLOTQLS 


ADJSMX 
ADVANC** 
AXLGLG** 
AXLI LI * * 


"CHKWRD 

EOF* 


CARDIN < 


IRDCRD+ 

NAMFIL+ 

L. 


EX ITPL* * 

FETADD+ 

FETDEL+ 

FORM** 

FTNBIN* 


* indicates a routine in the FORTRAN subroutine library. 

** indicates an NPS routine in the BCS FORTRAN subroutine 
library. 

*- indicates a routine in the DYLOFLEX alternate subroutine 
library. 

All others are local to L221 (TEV156). 



I 


Table 4. —(Continued) 

OVERLAY (L22 1,3,0) 

PROGRAM PLOTQLS 


GDLGLG** 

GDLI LI * * 

GMNMX 

INITIR+ 

LOCF* 

NOLGB** 

NOLGL** 

NOSLIB** 

NOSLLLI** 

PFLGLG* * 

PFLILI* * 

REDSOL ^READTP* 
REQFL+ 

SC4 20N* * 

SECOND* 

STARTR* 
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Table 4.— (Concluded) 

OVERLAY (L22 1,3,0) 

PROGRAM PLOTQLS 


STCHSZ * * 
STFONT* * 
STLNOR** 
STLNST** 
STNDIV* * 
STNCHR* * 
STNPTS* * 
STSIGF* * 
STSPEC** 
STSUBJ * * 
STS 20 B** 
STTOOL* * 
STTXTR** 


TITLEG** 



I/O Devices of PLOTQLS 

The I/O devices used in PLOTQLS are shown in figure 10. For a complete description of 
the input data cards and files, see sections 6.3 and 6.4 in volume I of this document. 



TAPE99 


Figure 10.— Overlay (L221 ,3,0) PLOTQLS Input/Output 
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3.5 DATA BASES 


The L221 (TEV156) data bases include I/O magnetic files (either tape or disk), internal 
scratch (temporary) storage files, and labeled common blocks. 

3.5.1 INPUT DATA 

The input data for L221 (TEV156) is in two forms; data cards and magnetic files (either 
tape or disk). 


Card Input Data 

For a complete description of the L221 (TEV156) card input data, see section 6.3 in 
volume I of this document. 


Magnetic Files (Tape or Disk) 

For a complete description of the L221 (TEV156) disk or tape input data, see section 6.4 
in volume I of this document. The possible files read are IFTAPE, INTAPE, and 
LDTAPE. 

3.5.2 OUTPUT DATA 

The output results of L221 (TEV156) may be of three types: printed results, magnetic 
files, and/or punched cards. 


Printed Output Data 

For a complete description of the printed output data, see section 6.5.1 in volume I of 
this document. 


Magnetic Files (Tape or Disk) 

For a complete description of the magnetic file output data, see section 6.5.2 in volume 1 
of this document. The possible output files are IPLTPE, IRTAPE, and NEWTPE. 

Punched Card Output Data 

For a complete description of the punched card output, see section 6.5.3 in volume I of 
this document. 

3.5.3 INTERNAL DATA 

L221 (TEV156) uses two methods to pass data between sections of the program: labelled 
common blocks and scratch (temporary) magnetic files (either tape or disk). 

Magnetic Files (Scratch Disk Files) 

L221 (TEV156) uses two disk files, SCRATCH and SCRAT2, for the temporary storage 
of data. Both files, SCRATCH and SCRAT2, will be returned to the system through the 
use of subroutine RETURNF 1 when L221 (TEV156) execution is terminated. 


Clemmons, R. E.: Programming Specifications for Modules of the Dynamic Loads System to 
Interface with FLEXSTAB. NASA Contract NAS1-13918, BCS-G0701, September 1975. 

(Internal Document.) 
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SCRATCH 


SCRATCH is a temporary scratch file written and read by the 1,0 primary overlay, 
FINDRMS. The file contains the power spectral density (PSD) solution frequencies and 
responses required to calculate loads. The data for each frequency is written as a single 
record with a standard FORTRAN binary WRITE statement. The contents of SCRATCH 
are displayed in figure 11. 


Array size 


NDOF x 1* 

l Q l i 


NDOF x 1 

IQMAGR) j 


r 

GFSj 


i* 

FFS, 


i* 

FYDAMP; 

j 


End-of-File 






Repeat for 
i = 1, NFREQ 


where: FFS 

FYDAMP 

GFS 

NDOF 

NFREQ 

Q 

QMAGR 


= The computed Kussner lift growth function for frequency i 
= The computed feedback transfer function for frequency i 
= The computed Wagner lift growth function for frequency i 
= The number of degrees of freedom 
= The number of frequencies 

= The array of generalized coordinates for frequency i 

= The array containing the magnitudes of |Q) 's elements for 
frequency i 

Signifies that arrays are complex and twice the size indicated 


Figure 1 1.— Contents of SCRATCH— A Temporary Scratch File 
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SCRAT2 

SCRAT2 is a temporary scratch file written and read by the 1,0 primary overlay, 
FINDRMS. The file contains the interpolated tabular input spectrums to be used in the 
calculation of {A} and {N 0 }. The arrays are written onto SCRAT2 with standard 
FORTRAN binary WRITE statements. The contents of SCRAT2 are displayed in 
figure 12. 


Array size 
NELE 


where: 


(SPFRlj 


End-of-File 



► 


Repeat for 
i = 1, NFREQ 


NELE = The number of unique spectrum. elements per frequency, 

(NFORC* (NFORC+1 ))/2 

NFREQ = The number of frequencies 

NFORC = The number of excitation forces 


Figure 12.— Contents of SCR AT2— A Temporary Scratch File 
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Common Blocks 


Table 5 displays the common blocks used in the program and the programs in which 
they are used. 

The labeled common blocks are used for communication between the main and primary 
overlays, and for communication between routines in a primary overlay. The block 
names and contents are described on the following pages. 

Blank common is used in the primary overlays of L221 (TEV156) as a variable length 
working storage area. Each overlay calculates the core required to process the problem 
being analyzed, and calls subroutines of VARDIM (variable dimensioning storage) 
subroutines to perform the bookkeeping of array allocation. Pointers to the array 
locations are stored in labeled common blocks. 


Table 5 .— Common Blocks Used in Each Overlay 
Common Blocks 



ALOVER 

BUFFRl 

BUFFR3 

BUFFR8 

BUFF11 

BUFR12 

CASLOC 

CHKPRT 

CNSTNT 

COMLOC 

CORSIZ 

GNCLOC 

INOUT 

KDEPEN 

LODLOC 

NPSCMN 

NPSCOM 

NPSZQ1 

NUMERR 

PLTBND 

PLTDTA 

RWBUFF 

SIZOPT 

SRTDTA 

TAPBUF 

TAPNUM 

VARDIM 

blank 

OVERLAY (L2 21, 0, 0) 







PROGRAM L22 lvc 

* 

* 

♦ * 

***** 



OVERLAY (L2 21, 1, 0) 







PROGRAM FINDRMS 

***** 

***** 

***** 

* 

* * 

* * * 

OVERLAY (L221, 2, 0) 







PROGRAM SORTQLS 

* 

* * 

* * 

* 

* * * 

* * 

OVERLAY (L221, 3,0) 







PROGRAM PLOTQLS 

* 

* * 

# * 

***** 

* * * 

* * 


The tables on the following pages display all of the common blocks within L221 
(TEV156). Table 5 indicates which blocks are used by each of the overlays. Tables 6 
through 9 show (per overlay) which routines require the common blocks. Finally, the 
contents of each labeled common block are described in table 10. 

Description of the Labeled Common Blocks 

Table 10 describes each variable contained in all of the labeled common blocks of L221 
(TEV156). The common blocks are ordered alphabetically, with each variable being 
described according to its location within the common blocks. 


Table 6.— Common Blocks Used in the Routines of L221vc 


Common Blocks 


J E-* 

CQ 


z 

z 

*-t 

05 

Q 

1 05 

H 

6 

Z 

o 

o 

05 

Z 

■ CL 

CO 

u 

u 


W 

CQ 

1 * 

05 

O 

if) 

CO 

CO 

Z 

Eh 


o 

z 

CL 

CL 

CL 

D 


! u 

o 

w 

z 

z 

z 

Z 

CL 


OVERLAY (L221, 0, 0) 
PROGRAM L221 vc 


CHKWRD 

IRDCRD 

PRGBEG 

PRGEND 
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Table 7 .— Common Blocks Used in the Routines of FINDRMS 

Common Blocks 



Global 3 

Local to (L221, 

1,0) - FINDRMS 

ALOVER 

CORSIZ 

INOUT 

NUMERR 

BUFFRl 

BUFFR3 

BUFFR8 

BUFF11 

BUFR12 

CASLOC 

CHKPRT 

CNSTNT 

COMLOC 

GNCLOC 

KDEPEN 

LODLOC 

RWBUFF 

SIZOPT 

TAPNUM 

VARDIM 

blank 

OVERLAY (L221, 1,0) 






PROGRAM FINDRMS 

* * * * 

***** 

* * * * 

* * * * 


CONDAT 

* * 

* 

* * * 

* 

* * 

CORREL 

* 





DELETR 

* 


* * 


★ * 

FEDBAK 






FRCALC 

* * 





GCMAT 

* * 

* 

* * * * 

* * * 

* * 

GRDPEN 






INITIR 

* * 


* 


* * 

LOCATR 

♦ 


* * 


* * 

LOMAT 

* * 

* 

* * * 

* * * * 

* * 

MATINT 




* 


MAT RED 

* * 



* 

* 

OCJTAP 

* * 





READTP 




* 


SAVEQ 

* * ! 





SINTRP 

♦ 

* 




SOLVEL 

* 





SOLVEQ 

* * 





SPCTRM 

* * 





STARTR 

* * 


* 


* * 

STATEL 






SUMATS 






WRTETP 

VTTTT'T-J 




* 



a Global means the block is used in more than one overlay 


Table 8.— Common Blocks Used in the Routines of SORTQLS 


Common Blocks 


Global 3 


Local 

to 

(L2 

to 
1— * 

2, 


to 


i 

Eh 

o 


< 

fa 

) 

X 


w 

w 

Eh 

(X 

0$ 

o 

fa 

Eh 

D 


M 


> 

eg 

D 

H 

O 4 

(J 


Q 

CQ 


g 

c 

o 


o 

as i 

« 

5; 

m 

Eh 

& 


2 


9 

o 

z 

D 

X 

O 


fa 

< 


< 


< 

u 

w 

£ 1 

U 

u 

fa 

CO 

Eh 


> 

XI 


- SORTQLS 


OVERLAY (L221, 2,0) 



PROGRAM SORTQLS 

* * * * 

* * * 

CHKWRD 



DELETR 

* 

* * 

INITIR 

* * 

* 

KARDIN 

* * 

t 

LOCATR 

* 

♦ * 

READTP 


* 

REDSOL 

* 


STARTR 

* * 

* 

WRTETP 


• 

XFER 



a Global r 

neans the block is used 

in more than one overlay 




VARDIM 


Table 10.— Description of the Labeled Common Blocks 


LABELLED COMMON NAME: ALOVER 


DESCRIPTION: Variables used between the primary overlays 


NO. 


VARIABLE 


TYPE 


DIM. 


ENG. 
[NOMENCLATURE 


DESCRIPTION 


MDOF 


MFREQ 


MS TATE 


MPLTPE 


MRTAPE 


NSETS 


ISET 


NG ROUP 


10 


Contains the number of 
degrees of freedom for 
SORTQLS and PLOTQLS . 

Contains the number of 
frequencies for use in 
SORTQLS and PLOTQLS . 

Indicates which static- 
elastic solution type 
will be used. 

Contains the name of 
the plot file for use 
in SORTQLS and PLOTQLS 

Contains the magnetic 
file name whose con- 
tents are the load 
transfer functions to 
be used in SORTQLS 
and PLOTQLS . 

Number of sets of 
load equations. 

Array containing the 
number of loads in 
each set. 

Contains the number of 
items to be contained 
on the magnetic file 
defined by MRTAPE for 
each frequency to be 
used by SORTQLS and 
PLOTQLS . 


a Variable types are as follows: C = complex, H = holler i th, 
O = octal, R = real 


integer, L = logical, 
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i ■ i 


mu 


III! I 


I II 


Mill I 




Table 10.— (Continued) 


LABELLED COMMON NAME: BUFFRl 

DESCRIPTION : Buffer Area for the Generalized Coordinate 

Scratch File, 



■ 

— 

~ 


_ _ 

NO. 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 







1 

LAI 

I 



Length of ARAY1 

2 

ARAY1 

R 

1000 


Buffer area used for 
the generalized coordi- 
nates (scratch) . 
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Table 10.— (Continued) 


LABELLED COMMON NAME: BUFFR3 

DESCRIPTION: Buffer Area for the Generalized Coordinates 

and Loads. 







NO. 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 







1 

2 

LA3 

ARAY3 

I 

R 

1000 


Length of A RAY 3 

Buffer area used for 
the Generalized 
Coordinates and Loads. 
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Table 10.— (Continued) 


LABELLED COMMON NAME: BUFFR8 

DESCRIPTION: Buffer Area used for the Input Matrices 








NO. 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 







1 

LA8 

I 

1 


Length of A RAY 8 

2 

ARAY8 

i 

| 

R, • 

540 


Buffer area used for 
the input matrices. 
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Table 10.— (Continued) 


LABELLED COMMON NAME: BUFFll 

DESCRIPTION: Buffer Area used for Plot Data 







NO. 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 

1 

2 

LAll 

ARA11 

I 

R 

1000 


Length of ARAll . 

Buffer area used for 
plot data. 
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Table 10.— (Continued) 
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Table 10.— (Continued) 


LABELLED COMMON NAME ; CASLOC _____ 

DESCRIPTION: Contains the Locations of the Case-dependent 

arrays in blank common (VARDIM storage). 


VARIABLE 


TYPE 


ENG. 

NOMENCLATURE 


DESCRIPTION 


LCSO 


LDAMP 


LFREQ 


LINSP 


LNDEL 


Pointer to {CSO}, 
special feedback co- 
efficient matrix. 

Pointer to ^CSl 1 , 
special feedback co- 
efficient matrix. 

Pointer to {CS2}, 
special feedback co- 
efficient matrix. 

Pointer to {DAMP}, the 
damping coefficient 
matrix . 

Pointer to {FLL } , the 
distance between grad- 
ual penetration panels 
Used to calculate the 
time log. 

Pointer to {FREQ}, the 
frequency solution 
matrix. 

Pointer to {INSP }, the 
tabular input spectrum 
matrix . 

Pointer to {NDEL }, the 
degrees of freedom 
deletion matrix. 



Table 10.— (Continued) 
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Table 10.— (Continued) 


LABELLED COMMON NAME 


CHKPRT 


DESCRIPTION: 


Problem Printout Option 


NO, 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 

1 

! 

i 

ICKPRT 

I 

i 

! 


Problem Printout 

Option 

=0 The input matrices 
are not printed. 

=1 The input matrices 
will be printed 
as read. 

=2 The input matrices 
plus the inter- 
mediate data will 
be printed for 
the first 3 fre- 
quencies . 

=5 Special program 
checkout feature. 

It prints 1 and 2 
above plus a record 
of all the VARDIM 
routines called. 
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Table 10.— (Continued) 


LABELLED COMMON NAME: 


CNSTNT 


DESCRIPTION: In P ut Constants for FINDRMS 


VARIABLE 


NOMENCLATURE 


DESCRIPTION 


ALPHAl 
BETAl 
GAMMA 1 


Wagner indicial lift 
growth function co- 
efficients . 




ALPHA2 
BETA2 
GAMMA 2 


OCOEF 


Kiissner indicial lift 
growth function co- 
efficients . 


The scalar constant 
for the equation of 
motions and loads 
equation forcing func- 
tion coefficient 
matrices . 

The velocity (true air 
speed) . 

Exponential Input 
Spectrum power 
constant. 



Table lO.-(Continued) 


LABELLED COMMON 

NAME:- 

CNSTNT 

(Continued) 


DESCRIPTION: 



















TYPE 


ENG. 


NO. 

VARIABLE 

DIM. 

NOMENCLATURE 

DESCRIPTION 







18 

OEX 

R 


V 

Exponential Input Spec- 
trum power constant. 

19 

TSCALE 

R 


L 

The turbulence scalar 
constant . 

20 

GPSCAL 

R 


c 

The gradual penetra- 
tion scalar constant. 

21 

TABFS 

R 



The tabular input 
spectrum frequency 
scalar . 

22 

TABS 

R 

1 



The tabular input 
spectrum power scalar 
constant . 

23 

FDSCAL 

R 



Contains the scalar 
constant used to scale 






the frequencies in the 
Frequency- Dependent 
Frequency Matrix 
(FREQM}. 

24 

GPSCAL 

R 


c 

Contains the gradual 
penetration scalar 






constant of (C^K 

25 

G 

i 

R 


G 

The feedback gain co- 
efficient of the feed- 






back transfer function/ 
FYDAMP . 

26 

RLIMIT 

R 


6 r 

Special feedback 





lim 

rudder limit. 
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Table 10.— (Continued) 
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Table 10.— (Continued) 


LABELLED COMMON NAME : COMLOC 

DESCRIPTION: Specifies the last array entry found by L0CATR. 


NO. 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 

1 

JTHENT 

I 



Specifies the catalog 
(KATLOG) entry of 
the last array located 
by routine LOCATR , 

( LOCATR is a VARDIM 
routine) . 
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Table 10.— (Continued) 


LABELLED COMMON NAME: CORSIZ 

DESCRIPTION: Core size 







NO. 

VARIABLE 

TYPE 

. j 

DIM. 

ENG. 

NOMENCLATURE 

.. 

DESCRIPTION 







1 

— 

MAXCOR 

1 

0 



Contains the maximum 
allowable core size. 
The default of 
200000 8 may be changed 
via card input. 
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Table 10.— (Continued) 


LABELLED COMMON NAME : GNCLOC 

DESCRIPTION: Locations of the Generalized Coordinate Arrays 










TYPE 


ENG . 


NO. 

VARIABLE 

DIM. 

NOMENCLATURE 

DESCRIPTION 







1 

LPHIT 

I 



Pointer to ['J ] , grad- 
ual penetration summa- 
tion matrix used to 
calculate {C 3 }. 

2 

LC2 

I 



Pointer to { C 2 } - 

3 

LC3 

I 



Pointer to {C 3 }. 

4 

LSUMM 

I 



Pointer to [SUMM], the 
equations of motion 
summation matrix. 

5 

LSUMC 

I 



Pointer to {SUMC }, the 
forcing function matrix, 

6 

LMl 

I 



Pointer to [M^ ] , an 
equations of motion 






matrix . 

7 

LM2 

I 



Pointer to [M 2 1, an 
equations of motion 






matrix . 

8 

LM3 

I 


i 

Pointer to [M 3 ], an 
equations of motion 
matrix. 

9 

LM4 

I 



Pointer to [M 4 ] , an 
equations of motion 






matrix. 

10 

LM5 

I 



Pointer to [ M 5 ] , an 
(equations of motion 
matrix. 
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Table 10.— (Continued) 


LABELLED COMMON NAME: 


INOUT 


DESCRIPTION: Input/Output File Names 


NO. 


VARIABLE 


TYPE 


DIM. 


ENG. 

NOMENCLATURE 


DESCRIPTION 


1 


INFIL 


I 


Name of the file used 
for card input data, 
(TAPES) . 


2 


IUTFIL 


I 


3 


IPFIL 


I 


Name of the file used 
for printed output, 

( TAP E 6 ) . 

Name of the file used 
for punched output, 
(TAPE7 ) . 




Ltli] 



Table 10.— (Continued) 


LABELLED COMMON NAME 


KDEPEN 


DESCRIPTION: 


Locations of the arrays needed for multiple 
Frequency- Dependent Matrix Problems. 


NO. 


VARIABLE 


TYPE 


DIM. 


ENG. 

NOMENCLATURE 


DESCRIPTION 


1 


0MEG1 


R 


Frequency at which 
input matrices are 


defined. 


2 

3 

4 

5 

6 


LM4 1 
LM51 
LM61 
LC3 1 
LPHITl 


I 

I 

I 

I 

I 


7 


LFREQM 


I 


Pointer to [M41 ] , 

[M 4 ] at 0MEG1 . 

Pointer to [ M51 ] , 

[ M 5 ] at 0MEG1 . 

Pointer to [M61] , 

[Ml at 0MEG1 . 

D 

Pointer to {C31}, 

{ C 3 } at OMEG1 . 

Pointer to [PHITl], 
[?] or [£] 
at 0MEG1 . 

Pointer to {FREQM}, 
the frequencies (in 
CPS) at which [ Mj , 

[ M ] , [ M ] , and {C- } 
or D [J] defined. 


8 


0J4EG2 


R 


9 


LM4 2 


I 


10 


LM52 


I 


Frequency at which 
input matrices are 
defined. 

Pointer to [M42], 
[M 4 ] at OMEG2. 

Pointer to [M52], 
[M 5 ] at OMEG2. 
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Table 10.— (Continued) 


LABELLED COMMON NAME: 


LODLOC 


DESCRIPTION : Locations of the Load Set-Dependent Arrays 


VARIABLE 

TYPE 

LLOAD 

I 

LLOADC 

I 

LNLDC 

I 

LCORSU 

I 

LCORLS 

i 

LQMAGL 

I 

LSPEC 

I 

LA 

I 

LPM 

I 

LPHITB 

I 

LC2B 

I 

LC3B 

I 



ENG. 

NOMENCLATURE 


DESCRIPTION 


Pointer to { LOAD > , a 
correlation array. 

Pointer to [ LOADC ] , a 
correlation array. 

Pointer to { NLDC } , 
number of loads to be 
correlated with each 
load . 

Pointer to [CORSUMl , 
a correlation array. 

Pointer to [CORLST], 
a correlation array. 

Pointer to {QMAGR} , 
the magnitude of { Q > • 

Pointer to {SPEC}, 
the output spectrum. 

Pointer to {A}, the 
RMS values. 

Pointer to { N o } , the 
zero crossings. 

Pointer to [?]# 
the summation matrix 
for gradual penetra- 
tion calculation of 
(C 3 } . 

Pointer to {C 2 }« 
Pointer to {C^}. 











Table 10.— (Continued) 


LABELLED COMMON NAME 


LODLOC 


(Continued) 


DESCRIPTION: 








NO. 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 







13 

LSUMMB 

I 



Pointer to [ SUMM] , 
the load summation 
matrix. 

14 

LSUMCB 

| 

I 



Pointer to { SUMC> , 
the load transfer 
function. 

15 

LPHIOM 

I 



Pointer to { PHI0ME } , 
the spectrums at 

16 

LMBl 

I 


i 

Pointer to [M^J, the 
load matrix. 

17 

LMB2 

I 


i 

Pointer to [Mjl / the 
load equation matrix. 

18 

LMB3 

I 



Pointer to [M-] , the 
load equation matrix. 

19 

LMB4 

I 



Pointer to [M^l, the 
load equation matrix. 

20 

LMB5 

I 



Pointer to [ Hi ] , the 
load equation matrix. 

21 

LMB6 

I 



Pointer to [Mg] , the 
load equation matrix. 

22 

LSBl 

I 



Pointer to [ S, ] , a 
load feedback matrix. 

23 

LSB2 

I 



Pointer to [Sj]/ a 
load feedback matrix. 

24 

LSB3 

•I 



Pointer to [S ^] , a 
load feedback 5 matrix. 
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Table TO.— (Continued) 



1 


Table 10.— (Continued) 


LABELLED COMMON NAME: . NPSCMN 


DESCRIPTION: NPS routines buffer area 


VARIABLE 


NOMENCLATURE 


DESCRIPTION 


Scratch storage area 
used by the NPS 
plotting subroutines . * 


Boeing Computer Services, Inc., "Numerical Plotting System - 
Users Manual"? BCS Document BCS-G0509? March, 1976. 

(Internal document) 




Table 10. —(Continued) 


NPSCOM 

LABELLED COMMON NAME: 


DESCRIPTION: NPS routines 



Boeing Computer Services, Inc - , "Numerical Plotting System - 
Users Manual"; BCS Document BCS-G0509; March, 1976. 

(Internal document) 



Table 10.— (Continued) 


LABELLED COMMON NAME: NPSZQ1 

DESCRIPTION: NPS routines 







NO. 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 







1 

DUM3 

R 

• 

3 

i 

i 

Scratch storage area 
used by the NPS plot 
subroutines. * 


* Boeing Computer Services, Inc.; "Numerical Plotting System - 
Users Manual"; BCS Document BCS-G0509; March, 1976. 

(Internal document) 
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Table 10.— (Continued) 


LABELLED COMMON NAME: NUMERR 


DESCRIPTION: Number of Diagnosed Errors 


i ;• 







VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 

l 

NERROR 

I 

1 


Contains the number 
of fatal errors diag- 
nosed during the 
execution of L221 
(TEV156) . 

i 
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LABELLED COMMON NAME: PLTBND 


DESCRIPTION: Plot Identification 


NO. VARIABLE 


TYPE 


ENG. 

DIM. NOMENCLATURE 


DESCRIPTION 


1 I PLOT 


2 I DP LOT 


3 SPECS 


4 BUFFER 


Indicates the status 
of the plots. 

=0 plot has not been 
previously called. 
=1 plotter has been 
previously called . 
=2 close the plotter. 

Array containing the 
user's plotting in- 
formation (name, phone 
no. , mail stop, and 
organization) . 

Array containing the 
internal NPS plotting 
specifications . 

Buffer area for NPS 
plotting routines . 
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Table 10.— (Continued) 



0. VARIABLE 


ENG. 

TYPE DIM. NOMENCLATURE 


DESCRIPTION 


ITYPE 


IPAIR 


JCTATE 




IRRAY 


SCALX 


SCALY 


Contains the load set 
from which to plot. 

Option indicating the 
variable plotted. 
l=Q's 
2=Loads 
3=Spectrums 
4=Root Mean Square 
of Loads. 

Denotes whether the 
standard and the 
static-elastic solu- 
tions will be paired 
for plotting. 

=0 separate 
=1 pair 

Index on the solution 
to be plotted. 

=0 standard 
=1 static-elastic 

Number of responses 
to plot. 

Array containing the 
numbers of vectors 
to be plotted. 

Scalar constant for 
the frequency array. 

Scalar constant for 
the output responses. 




Table 
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l— (Continued) 


(Concluded) 


ENG. 

NOMENCLATURE DESCRIPTION 

Misc. title blocks 
for labeling the plot 
frames . 

Grid limit input 
=0 set to data limit 
=1 input on cards. 

Contains the X and Y 
grid limits 
GRIDL ( 1 ) = X min.lmts. 
GRIDL ( 2 ) =X max.lmts. 
GRIDL ( 3 ) =Y min.lmts. 
GRIDL ( 4 ) =Y max.lmts. 

Axis scale type 
=0 Linear-linear 
-1 log-log 

Array containing the 
X-axis plotting 
label . 

Array containing the 
Y-axis plotting 
label. 




Table 10.— (Continued) 


LABELLED COMMON NAME: RWBUFF 


DESCRIPTION: 


Buffer area for Routines READTP/WRTETP 




i 




NO. 

l 

VARIABLE 

.. 

TYPE 

DIM. 

ENG . 

NOMENCLATURE 

DESCRIPTION 


i — 





1 

ITELL 

H 



Contains 8HBUFFSI ZE 

2 

, IS I ZE 

I 



Length of DECRS . 

3 

i 

i 

i 

i 

i 

DECRS 

1 

i 

i 

1 

1 

10000 

1 

1 

1 

1 

1 

1 

i 

Buffer area used by 
1 the routines READTP/ 
WRTETP . 
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Table 10.— (Continued) 


LABELLED COMMON NAME: SIZOPT 

DESCRIPTION* Contains the problems , dimensions, options, and 
tape numbers to be used in FINDRMS. 







NO. 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 

1 

2 

3 

4 

5 

FREQL 

ICORR 

IDAMP 

IFDBAK 

IFREQ 

R 

I 

I 

I 

I 

10 


Frequencies at which 
LOAD/DOF will be 
printed . 

Number of loads to be 
correlated with others. 

Option which indicates 
the number of different 
structural damping 
factors . 

=0 the same factor 
for all freedoms. 

=1 a different damp- 
ing factor for 
each frequency. 

Option which chooses 
the feedback loop to 
be used. 

=0 no feedback loop 
used . 

=1 standard feedback 
loop used. 

=2 special feedback 
loop used. 

Option which chooses 
the frequency matrix 
input . 

=-l Frequencies will 

be read from cards 
in cycles/second. 
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Table 10.— (Continued) 


LABELLED COMMON NAME: SIZ0PT (Continued) 

DESCRIPTION: 







NO. 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 

6 

7 

8 

IFTAPE 

INTAPE 

INTZRO 

H 

H 

I 



=0 Frequencies will 
be read from cards 
in radians per sec 

=1,2,3, 4, 5 The fre- 
quency matrix will 
be generated 
according to 
Appendix A. 

Denotes the location 
of the Tabular Input 
Spectrums. They will 
be read from cards if 
blank or will be read 
from tape "IFTAPE" if 
defined . 

Denotes the location 
of Input Matrices. 

They will be read from 
cards if blank or will 
be read from tape 
"INTAPE" if defined. 

Option which chooses 
the starting point 
for the RMS integra- 
tion. 

=0 start at zero. 

=1 start at FREQ^ . 
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Table 10.— (Continued) 


LABELLED COMMON 

NAME: 

SIZOPT 

(Continued) 


DESCRIPTION: 













■MM 




n 


TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 




■MMH 



9 

10 

11 

IPLTPE 

IPLRMS 

IPRINT 

i — 

H 

I 

I 

1 


Indicates whether or 
not a plotting file 
will be created. Plott- 
ing information will 
not be written onto 
a file if IPLTPE is 
blank, otherwise it 
will be written onto 
tape "IPLTPE" if de- 
fined . 

Option requesting that 
the 5 values be 
written on the plot 
tape. 

=0 A values not saved. 
=1 A values saved on 
"IPLTPE" for each 
frequency . 

Option which chooses 
the type of printing 
to take place. 

=0 print only the 

final A values and 
the number of zero 
crossings. 

=1 print the general- 
ized coordinates 
and load transfer 
functions for each 
frequency as well 
as "0" above. 
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Table 10.— (Continued) 


LABELLED COMMON 

NAME : _ 

SIZOPT 

(Continued) 


DESCRIPTION: 













NO. 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 







12 

IPUNCH 

I 



Option which chooses 
whether or not the 
final a values and 
the number of zero 
crossings will be 
punched on cards. 

=0 no punched output 






=1 punched output. 

13 

IREDUC 

I 



Not used. 

14 

IRMSPR 

I 



Option which chooses 
the type of printing 
for the A values. 

=0 print only the 
final A values. 






>0 print the "A 

value for each of 
the last IRMSPR 
frequencies . 

15 

IRTAPE 

H 



Denotes whether or not 
the Q 1 s and the load 
transfer functions 
will be written. The 
information will not 
be saved if left blank, 
otherwise the infor- 
mation will be written 
onto tape "IRTAPE" . 

16 

I-SPEC 

I 

I 


Option which chooses 
the type of input 
spectrum to be used. 
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Table 10.— (Continued) 


LABELLED COMMON NAME: 


SIZOPT (Continued) 


DESCRIPTION: 


NO, 


VARIABLE 


TYPE 


DIM. 


ENG. 

NOMENCLATURE 


= 0 
= 1 
= 2 


DESCRIPTION 


Dryden 
Von Karman 
Expoential 
NSP Tabular input 
spectrum defined 
at the NSP point 


17 


ISTATE 


I 


Option which chooses 
the static-elastic 
solution . 

=0 no static-elastic 
solution. 


18 


LABEL 


H 


7 


=1 static-elastic 
solution found 
with rigid body 
degrees of free- 
dom in col.l & 2. 

=2 static-elastic 
solution found 
after zeroing 
specified columns . 

Contains miscellaneous 
labeling information 
for labeling the data 
case. 


19 


LDTAPE 


H 


Denotes the location 
of the load equations. 
If the load equations 
and/or f requency-de- 
pendent load equations 
are contained on tape 
11 INTAPE" or if they 



Table 10.— (Continued) 


LABELLED COMMON NAME: SIZ0PT (Continued) 

DESCRIPTION: 







NO. 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 







20 

21 

22 

23 
2 4 

LPDOF 

NDOF 

NDOFD 

NFORC 

NFREQ 

I 

I 

I 

I 

I 



are to be read from 
cards the variable 
LDTAPE is to be left 
blank . 

Option which indicates 
whether or not the 
loads per degrees of 
freedom will be found. 
=0 will not be found. 

>0 loads due to each 
degree of freedom 
will be found and 
printed. 

Contains the number of 
degrees of freedom. 

Contains the number of 
degrees of freedom to 
be deleted. 

Contains the number of 
excitation forces for 
the handling of multi- 
ple forcing functions. 

Contains the number of 
frequencies over which 
the solution will be 
found . 
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». -(Continued) 




Table 10.— (Continued) 


LABELLED COMMON NAME: SIZ OPT ( Co nclud ed) 

DESCRIPTION: 


;o. 


VARIABLE 


TYPE 


DIM. 


ENG. 
[NOMENCLATURE 


DESCRIPTION 


Nulma t 


Matrix to be 
Read 



EOM 


Loads 

1 

[M 1 ] 

or 

[Mj| 

2 

[m 2 ] 

or 

[m 2 ] 

3 

[m 3 1 

or 

[m 3 ] 

A 


or 

[m a ] 

5 

[m 5 ] 

or 

[m 5 ] 

6 

[m 6 3 

or 

[m 6 ] 

7 

{C 2 } 

or 

{c 2 ) 

8 

{C } 

or 

{c 3 } 

9 




10 

[0] 

or 

M 

11 

[Sil 

or 

ISy] 

12 

[ S 2 J 

or 

[S 2 ] 

13 

[S 3 ] 

or 

[ s 3 ] 

14 

t s 4 ] 

or 

[S 4 ] 

15 

[S 5 ] 

or 

[Sg] 

16 

l S 6l 

or 

IS 6 ] 


Table W.-(Continued) 




ENG. 

NOMENCLATURE 


DESCRIPTION 


JSTATE 


IRRAY 


SCALX 


SCALY 


Index of the set of 
loads to be sorted. 

Option indicating the 
variable to be sorted. 

1 = Q 1 s 

2 = Loads 

3 = Spectrums 

4 = RMS load values 

Denotes whether the 
standard and the 
static- elastic solu- 
tions will be paired. 
=0 separate 
=1 pair 

Index on the solution 
to be sorted. 

=0 standard 
=1 static-elastic 

Number of responses 
to sort . 

Array containing the 
response numbers to 
be sorted. 

Scalar constant for 
the output fre- 
quencies . 

Scalar constant for 
the output responses. 
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Table 10.— (Continued) 
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Table 10.— (Continued) 


LABELLED COMMON NAME: TAPNUM 

DESCRIPTION: File Constants 








NO. 

VARIABLE 

TYPE 

DIM* 

ENG. 

NOMENCLATURE 

DESCRIPTION 

1 

2 

3 

— — — 

ITAPE 

IPR 

IRLCMP 

H 

I 

I 



Contains the magnetic 
file name (either 
” IN TAPE" or "LDTAPE" , 
if defined) from which 
MATRED is to read a 
matrix. 

Variable which con- 
tains the special 
trouble-shooting print- 
out option. 

Indicates the matrix 
type. 

=0 Real matrix 
=1 Complex matrix 
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Table lO.-(Continued) 



VARIABLE 


ENG, 

NOMENCLATURE 


DESCRIPTION 


NENTRY 


LWAVAL 


LWUSED 


MAXUSD 



Maximum number of 
arrays possible in 
VARDIM storage. 

Number of arrays 
currently defined in 
VARDIM storage. 

Last word available 
in blank common. 

Last word currently 
being used in blank 
common by a VARDIM 
array. 

First word address 
of KATLOG in blank 
common. 

Maximum core used 
by VARDIM since 
calling routine 
STARTR. 




Table 10.— (Concluded) 


LABELLED COMMON NAME: blank 

DESCRIPTION: Blank COmmon 







NO. 

VARIABLE 

TYPE 

DIM. 

ENG. 

NOMENCLATURE 

DESCRIPTION 

. 






1 

D 

R 

1000 

1 

j 


The blank common will 
be used to store the 
matrices under a vari- 
able dimensioning 
storage scheme 
(VARDIM) . 
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Description of Blank Common Used 

The array or matrix sizes required in the modules of L221 (TEV156) vary widely from 
one problem to another. In fact, in most cases more than 50% of the matrices are null 
and need not be stored. For these reasons the VARDIM routines described later in this 
section are used to dynamically allocate core, compact storage, and keep track of most 
array locations and sizes. 

The arrays independent of problem size are still defined in the standard FORTRAN 
fashion. Only the arrays using VARDIM will be discussed in this section. 

Program L221vc 

No VARDIM array storage is required. 

Program FINDRMS 

Three different blocks of arrays are defined in VARDIM storage as FINDRMS is 
executed. The arrays will be defined only if the statement under the heading 
option-dependent is true. 

• Block I is defined in the subroutine CONDAT and contains the arrays of input 
constants and frequencies (table 11). Block I is made up of card sets 3 through 11. 

• Block II is defined in the subroutine GCMAT and contains the arrays necessary to 
calculate the generalized coordinates (table 12). Block II is made up of card sets 12 
through 30. 

• Block III is defined in the subroutine LOMAT and contains the arrays necessary to 
calculate the load frequency response and the RMS values (table 13). All arrays of 
Block II except Q, QMAGR, and FLL are deleted before Block III is defined. 
Block III is made up of card sets 31 through 49. 
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Table 1 7.— Input Constants and Frequencies of Block / 


Array Name 

Size 


Option- dependent 

NSTEL 

NDOF 


I STATE>0 

NDEL 

NDOFD 


NDO FD> 1 

DAMP 

NDOF 


— 

CS2 

NDOF 



CS 1 

NDOF 

► 

I FDBAK=2 

CSO 

NDOF 



FREQ 

NFREQ 


-- 

FREQS 

NS P 



SPECIN 

NSP* 2*NELE 

I S P EC < 0 

S P FR 

2*NELE 


ISPEC<0 


where NELE 

= (NFORC* (NFORC <■ 1 ) ) /2 


NSP 

= ABS(ISPEC) 
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Table 12.— Arrays Contained in Block II 


Array Name 

Size 


Option- dependent 

Q 

2 * N DO F* N FO RC 

-- 

QMAGR 

NDOF+NFORC 

-- 

FREQM 

NKVAL 

NKVAL>1 

FLL 

NPAN 


NPAN>0 

PH IT 

2*NDOF*NPAN 

NPAN>0 

C2 

NDOF*NFORC 

-- 

C3 

2*NDOF*NFORC 

-- 

SLJMM 

2*NDOF*NDOF 

-- 

SUMC 

2*NDOF*NFOC 

— 

Ml 

NDOF*NDOF 

NULMAT (1 ) >0 

M2 



NULMAT (2) >0 

M3 



NULMAT (3) >0 

MU 



NULMAT (4) >0 

M5 



NULMAT (5) >0 

M6 

NDOF*NDOF 

NULMAT (6) >0 



IFDBAK>0, and 

SI 

NDOF 3 

►NDOF 

NULMAT (7) >0 

S2 



NULMAT (8) >0 

S3 



NULMAT (9) >0 

S4 



NULMAT (10) >0 

S 5 


f 

NULMAT (1 1 ) >0 

S6 

NDOF*NDOF 

NULMAT (1 2) >0 


Table 12.— (Concluded) 


Arr ay Name 


M4 1 

M42 

M51 

M52 

M61 

M62 

C31 

C32 

PHIT1 

PHIT2 


} 

} 

} 

} 


Size 


Opt ion- dependent 


NK VAL> 1 9 and 


NDOF*NDOF 


NDOF*NDOF 
2 *NDOF*NFORC 
2*NDOF*NFORC 
2*NDOF*NPAN 
2 *NDOF*NPAN 


{ 

{ 

{ 

{ 

{ 


NULMAT (4) >0 


NULMAT (5) >0 


NULMAT (6) >0 


NPAN <0 


NPAN>0 



Table 13.— Arrays Contained in Block III 


AlEay_Name 

LOAD 

LOADC 

NLDC 

CORSUM 

CORLST 

QMAGL 

SPEC 

A 

PM 

C2B 

PHITB 

C3B 

SUMMB 

SUMCB 

PHIOME 

MB1 

MB2 

MB3 

MB 4 

MB 5 

MB6 


Si ze 
ICORR 
ICORR+1 0 
ICORR 
ICORR* 1 0 
ICORR*1 0 
NLD*NFORC 


NLD*NFO RC 
2*NLD*NPAN 
2 *NLD*N FORC 
2*NLD*NDOF 
2 *NLD 

(NFORC* (N FORC + 1 
NLD*NDOF 


t 

NLD*NDOF 


QEtiEDzde nendent 
ICORR>0 


NPAN>0 


) / 2 -- 

NULMAT (1) >0 
NULMAT (2) >0 
NULMAT (3 ) >0 
NULMAT (4) >0 
NULMAT (5) >0 
NULMAT (6 ) >0 



Array N am e 


SB1 
SB2 
SB3 
SB4 
SB5 
SB 6 


M4 1 

M42 

M51 

M52 

M61 

M62 

PHIT1 

PHIT2 

C31 

C32 


} 

} 


Table 13.— (Concluded) 
Size 


NLD*NDOF 


NLD*N DO F 


NLD*NDOF 


NLD*NDOF 
2 *NLD*NPAN 
2*NLD*NPAN 
2*NLD*NF0RC 
2*NLD*NFORC 


Opti on -dependent 
IFDBAK>0, and 
NULMAT (7) >0 
NULMAT (8) >0 
NULMAT (9 ) >0 
NULMAT (10) >0 
NULMAT (1 1) >0 
NULMAT (1 2) >0 
NKVAL> 1 , and 

^ NULMAT (4) >0 

' NULMAT (5) >0 

i 

f NULMAT (6) >0 


NPAN>0 


| NPAN<0 
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Program SORTQLS 

Two or three arrays will be defined in VARDIM for SORTQLS. 

Array name Size Option-dependent 


FREQ NFREQ 

SOLI (NFREQ* NR) 

SOL2 (NFREQ*NR) IPAIR>0 


where NR is NDOF or NLD, depending upon the response being plotted (NDOF for 
QMAGR and NLD for QMAGL and/or SPEC), and IPAIR is 1 if pairing of standard and 
static elastic solutions has been requested. 

Program PLOTQLS 

The required VARDIM storage is the same as described in program SORTQLS above. 

Variable Dimensioning Storage Scheme (VARDIM) 


Analyses requiring the storage of matrices can lead to inefficient large core programs 
(i.e., coded for the maximum matrix sizes) unless some method is used to make the 
storage required dependent upon the individual problem’s size. A series of five 
subroutines (STARTR, INITIR, LOCATR, DELETR, and XFER) collectively known as 
VARDIM, has been written to handle the allocation of matrix storage during program 
execution. VARDIM uses blank common for all array storage. This method is possible 
on operating systems that place blank common at the end of all other program storage. 
The user must request enough field length (core) to provide sufficient array storage 
between the beginning of blank common and the end of the declared field length. 

VARDIM is not a matrix language. It does not provide matrix operations, only array 
storage allocations and the necessary bookkeeping. The subroutines were designed to be 
called by FORTRAN programs that do their own matrix storage (input/output), perform 
their own matrix operations, and calculate their own array subscripts. 

Some general features of the VARDIM scheme follow: 


1. Array storage may be variable within a single program run as well as between 
runs. Arrays may be defined or deleted at any time. It is not necessary to define 
them all at the beginning of the program. 

2. Each array will be identified by a six-character hollerith name (left-justified), may 
have from one to three dimensions, and may be a real or integer variable type. 

3. Newly defined arrays will be null. 
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4. The array storage is always compacted to use the first words of blank common. A 
newly defined array is always located after pre-existing arrays. If array i is 
deleted, then arrays i + 1 through n are moved forward to positions i through n-1. 

5. When any VARDIM routine is called, it checks to see that there are no duplicate 
array names or illegal dimensions and that core and/or the array catalog has not 
been exceeded. 

The bookkeeping performed by VARDIM is stored in a catalog array KATLOG and 
the labeled common block /VARDIM/. 


6. The VARDIM array catalog, KATLOG, is itself an array in VARDIM storage with 
the dimensions (6, NMAX, 1). Each array will have a six-word entry in KATLOG 
(i.e., one column). The six words contain: 

a. Name-6 hollerith characters, left-justified and blank filled 

b. Location-the first word address of the array in blank common 

c. Type-0 for integer and 1 for real 

d. Row dimension size 

e. Column dimension size 

f. Level dimension size 


7. All VARDIM routines except XFER contain the labeled common block /VARDIM/, 
which has the bookkeeping variables: 


NMAX 

NENTRY 

LWAVAL 

LWUSED 

LKAT 

MAXUSD 


Maximum number of arrays which may be defined for this program 
Number of arrays currently defined 

Last word available in blank common (length of blank common) 

Last word currently in use by VARDIM in blank common 

The first word address in blank common of the array named KATLOG 

Maximum core used by VARDIM since calling routine STARTR 


8. All VARDIM arrays except XFER contain the blank common definition 


COMMON D(l) 

DIMENSION ID(1), KATLOG (6,1) 
EQUIVALENCE (D,ID, KATLOG) 
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9. A program checkout feature is available. The VARDIM routines contain the 
following labeled common block. 

COMMON / CHKPRT / ICKPRT 

If ICKPRTssS, each VARDIM routine except XFER will print a message giving 
information about the array and core locations it is manipulating. There are two 
print lines for each message. 

NOTES: 

• For a description of the VARDIM routines used by the overlays of L221 (TEV156), 
see the discussion of the VARDIM routines beginning on the next page. 

• It will be necessary to use the REDUCE(-) control card to prevent the loader from 
reducing the program’s executable field length to dimensioned blank common 
size of 1. 


Discussion of VARDIM Routines 


The purpose and usage of each of the VARDIM routines will be briefly described. The 
variables below appear repeatedly in their argument lists. 


• NAMEA 

• LA 

• MA 

• NA 

• INTRA 

• LOCA 


An array name of up to 6 hollerith characters (left-justified) 

First dimension of array NAMEA 
Second dimension of array NAMEA 
Third dimension of array NAMEA 

A key to the VARDIM routines indicating NAMEA’s variable type 
(0 for integer and 1 for real) 

The location of array NAMEA-the first word address in blank 
common 


Subroutine STARTR 

STARTR must be called before any other VARDIM routine. It creates the array catalog, 
sets the maximum number of arrays, and determines the maximum array storage 
available. 

Calling sequence of STARTR: 

CALL STARTR (IMAX, LENGTH) 

Input 

IMAX The maximum number of arrays to be defined at any one time; 

if IMAX - 0, STARTR will assume 100 

LENGTH The number of words available in blank common for use by the 

VARDIM routines; if length is input as zero, STARTR will find 
the length available (difference between the beginning of blank 
common and end of field length) 
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Output 

LENGTH The number of words available in blank common; if LENGTH 

is returned as zero, there was not enough blank common to 
create the catalog array, KATLOG 

Subroutine INITIR 

INIT1R is called to initialize or define a new array. It must be called before any other 
VARDIM routine refers to the array. INITIR allocates the array storage, zeros the area, 
and makes an entry in the array catalog. When the array being defined already exists, 
the elements are simply set equal to zero if the dimensions are to remain the same. If 
the array size is to be changed, the old array is deleted and a new one defined. 

Calling sequence of INITIR: 

CALL INITIR (NAMEA,LA,MA,NA,INTRA,LOCA,IRI) 

Input 

NAMEA 
LA 
MA 
NA 

INTRA 
Output 

LOCA See subroutine argument definitions 

IRI Error code 

= 0, no error detected 

= -1, previous array of same name deleted before defining a 

new one 

= 1, maximum number of catalog entries was exceeded 

= 2, one of the arrays dimensions is 0 

= 3, blank common storage exceeded 

= 4, routine STARTR was not called beforehand. 

Subroutine DELETR 

DELETR is called to eliminate an array from VARDIM storage. The array’s entry in the 
VARDIM catalog will disappear and both the catalog and array storage compacted. 

Calling sequence of DELETR: 

CALL DELETR (NAMEA, IRD) 

Input 

NAMEA See subroutine argument definitions 


> See subroutine argument definitions in discussion 
I of VARDIM routines 
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Output 

IRD 


Error code 


= 0, no error detected 

= -2, NAMEA is not in the VARDIM catalog 


Subroutine LOCATR 


LOCATR is called to determine an array’s size, type, and location. LOCATR should be 
called just before handling the array because the array’s location changes as other 
arrays are deleted and added to VARDIM storage. 


Calling sequence of LOCATR: 

CALL LOCATR (NAMEA,LA,MA,NA,INTRA,LOCA,IRL) 
Input 

NAMEA See subroutine argument definitions 


Output 

LA 

MA 

NA 

INTRA 

LOCA 

IRL 


\ 


See subroutine argument definitions 


Error code 

= 0, no error detected 

= -1, no entries in catalog 

= -2, NAMEA is not in the catalog 


Subroutine XFER 


XFER is a subroutine called by DELETR to quickly move blocks of core. 
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4.0 EXTENT OF CHECKOUT 


Nine different data cases were used to verify that L221 (TEV156) runs correctly. The 
results of the nine cases were compared against answers produced by hand calculations 
and previous versions of the program (see footnote 1 in sec. 2). 

Checkout Problems 


No. Type of problem 

1 Standard execution with input equations of motion and load equations from cards 

2 Standard execution with equations input on a magnetic file; SORTQLS and 
PLOTQLS also 

3 Execution with frequency-dependent equations of motion and load equations input 
on cards 

4 Same as 3 except equations input via a magnetic file 

5 Multiple forcing functions input on cards 

6 Same as 5 except equations input via a magnetic file 

7 Same as 1 except freedoms physically deleted from cards; serves as a check case for 
number 8 

8 Same as 1 except freedoms deleted through use of a program option 

9 Magnetic file input of matrices (same as 4) with DYLOFLEX header arrays on the 
file 

The options used in the checkout problems are displayed in table 14. 


Boeing Commercial Airplane Company 
P.O. Box 3707 

Seattle, Washington 98124 
May 1977 
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Table 14.— Options Used in Checkout Cases 



OPTIONS 


NDOF 

NFREQ 

NPAN 

NKVAL 

NFORC 

NDOFD 

IFREQ 

ISPEC 

I STATE 

IFDBAK 

I DAMP 

INTZRO 

IPRINT 

IRMSPR 

IPUNCH 

ICKPRT 

IPLRMS 

LPDOF 

INTAPE 

IPLTPE 

IRTAPE 

IFTAPE 

Z 

VEL 

L(TSCALE) 
C(GPSCAL) 
TABFS 
TABS 


7 7; 7 7 

10 10?10 10 

0 5 ; 0 0 

0 0 ; 0 10 

0 0 ; 0 0 

0 0; 0 0 

0 0 ; 0 0 

1 1 ? -10 1 

0 2 ; 0 0 

0 1 ; 0 0 

0 1 ; 0 0 

1 0; 1 1 

1 1; 0 0 

3; 0 0 

0; 0 0 

2; 0 0 

1; 0 0 

2 ; 0 0 

INTAPE; - 

ONEll TWOll ; - THREll 
TAPE 3 ; - 


12 . 12 • ; 12 . 12 . 

2613. 812613. 83;2613. 83 
2613.83 


2500. 2500; 250D 2500. 


3 

4 

5 

6 

7 

8 

2; 7 

2; 7 

10 

210 

3 ; 72 

3 ; 7 2 

0 

0 

o 

O 

0; 0 

10 

20 

0; 0 

0; 0 

0 

0 

2; 6 

2; 6 

0 

0 

0; 0 

o 

o 

0 

0 

-l; o 

-1;0 

1 

1 

- 3 ? - 2 6 

- 3 ; - 2 ( 


1 . ? 0 . 
0 . ; 0 . 


0 

0 

0 

1 

0 

0 

0 

2 

0 

0 

EOMTPE 

IRTAPE 

1 . 

829.53 

2500. 

0. 

0 . 




1 . ; 1 . 
1 . ; 1 . 

1 . ; 1 . 
0 . ; 0 . 
0 . ; 0 . 
0 . ? 0 . 


8 

210 

0 

20 

0 

0 

0 

1 

0 

0 

0 

1 

0 

0 

0 

2 

0 

0 

DYFEOM 


; IFTAPE 

1 . ; 1 . 12 . 12 . 1 . 

1 . ; 1 . 2613 . 8 3 2613.82 829.53 


1- Jl- 2500. 2500. 2500. 



Table 14.— (Concluded) 
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16 Abstract 

L221 (TEV156) is a digital computer program available for execution on the CDC 6600. 
The program is capable of calculating steady-state solutions for linear second-order 
differential equations due to sinusoidal forcing functions. From this steady-state 
solution, generalized coordinates and load frequency responses may be determined. 
Statistical characteristics of loads for the forcing function spectral shape may also be 
calculated using random harmonic analysis techniques. 

The particular field of application of the program is the analysis of airplane response 
and loads due to continuous random air turbulence. Optional capabilities include 
frequency dependent input matrices, feedback damping, gradual gust penetration, 
multiple excitation forcing functions, and a static elastic solution. 

Program usage and a brief description of the analysis used are presented in volume I of 
this document. Volume II contains a description of the design and structure of the 
program to aid those persons who will maintain and/or modify the program in the 
future. 
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